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Abstract 

D ' We have theoretically studied the collective modes of a quasi two-dimensional (Q2D) Bose con- 

densate in the large gas parameter regime by using a formalism which treats the interaction energy 
I beyond the mean-field approximation. In the calculation we use the perturbative expansion for 

the interaction energy by incorporating the Lee, Huang and Yang ( LHY) correction term. The 
■ results show that incorporation of this higher order term leads to detectable modifications in the 

o , 

^ O^ , mode frequencies. 

^ , PACS numbers: 

in 

03.75.Kk Dynamic properties of condensates; collective and hydrodynamic excitations, superfluid flow. 

o ■ 

g . 03.75.Hh Static properties of condensates; thermodynamical, statistical, and structural properties. 

o 



o 



X 



05.30.Jp Boson systems (for static and dynamic properties of Bose-Einstein condensates). 
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I. INTRODUCTION 



Bose-Einstein Condensation (BEC) of dilute atomic gases has been achieved in a variety 
of atomic samples using different magnetic and optical traps 

mm 

. These advances 

have resulted in a investigations of properties of these ultra-cold gases in different conditions. 
At present, there is considerable interest in BEC in low dimensions where particle motion 
is restricted in one or more dimensions 

Q, Q, Q- This is because a 
Bose gas confined in low dimensions shows remarkably different behaviour than it shows in 
three dimensions. For example, a two-dimensional homogeneous Bose gas does not undergo 
Bose-Einstein condensation on cooling. Instead, there is a phase transition called Berezinski- 
Kosterlitz-Thouless transition in which system becomes superfluid without having long-range 
ordering. For a homogeneous Bose gas confined in one dimension, BEC does not exist even 
at absolute zero. However, in inhomogeneous confining potentials BEC can exist in lower 
dimensions also [isj ]. 

Recently, several groups have experimentally investigated the properties of Bose conden- 
sates confined in the low dimensions 0, S|- To investigate the macroscopic quantum 
behaviour of these systems, the study of collective modes has proved to be a fundamental 
tool. Usually, the theoretical description of these dynamical properties such as collective 
modes of a Bose condensate can be accomplished by Gross-Pitaevskii (CP) equation taking 
interaction energy within the mean-field approximation. However, it is known that the mean- 
field theory gives accurate results only in the regime of small values of gas parameter (pa^, p 
is density and is scattering length). With the recent experimental advances in this field, 
access to the large gas parameter regime has become feasible either by increasing the number 
of atoms in the condensate or by increasing the scattering length via Feshbach resonances. 
Therefore it is essential to assess the CP theory in this regime. In the present theoretical 
studies, we have calculated the collective mode frequencies of a Q2D gas in the regime of 
large gas parameter by using a perturbative expansion for the interaction energy. We have 
incorporated the Lee, Huang and Yang ( LHY) correction term in the energy expansion 
to go beyond the mean- field approximation |l6|. We find that detectable changes in the 
collective mode frequencies can occur when we incorpoarte the LHY correction term in the 
energy expansion. 



2 



II. GROSS-PITAEVSKII THEORY 



The GP equation for a condensate admits formulations in terms of the Hamilton least 
action principle with the action functional given as, 

S = I Ldt and L = [ Cdr. (1) 



Here C is Lagrangian density and is given by, 



where Vext is trapping potential, m is mass of the atom, e{p) is self interaction energy, ip is 
condensate wave function and p is densityd?/'!^). The normalization condition is, 

j\i.\^ir = jpir = N (total number of atom.). (3) 

By evaluating Lagrangian for a trial wave function, the solution of Lagrange's equations 
gives the parameters used in the trial wavefunction. 

In the large gas parameter regime, we use the perturbative expansion Q| for the in- 
teraction energy e{p) beyond mean field approximation by including the correction term 
calculated by Lee, Huang and Yang ( LHY) [l7| as following, 

e{p) = ^-^p[l + C,{palf% (4) 
m 

where Ci = 128/(15y^). Here the first term in the expression of interaction energy e{p) 
represents the mean-field energy whereas the second term represents the LHY correction 
term. It is also known that this expansion is valid upto pa^ ~ 10^^ 3, 18|. 

Further, we assume that the Bose condensate is trapped in a harmonic potential of the 
form, 

Vext = ^miuy + jy + ujIz""), (5) 

where ojx-.ojy.ujz are the trap frequencies along x, y and z-axes respectively. 

In the Q2D regime of confinement, the condensate is supposed to be confined in a trap 
which has frequency along one axis (z-axis) much larger than the frequencies along the other 
two axes (x-axis and y-axis). This implies that uj^ » ujx,ujy and motion in z-direction is 
frozen. However, the oscillator lengths aj = (/l/mtUj)^/^, (i = x,y,z) along the three axes 



of the trap are still much larger than the scattering length so that the collisions are three 
dimensional. For the scattering length larger or comparable to the oscillator length along 
the confining axis, the scatteri ng prop erties become significantly different from those of a 



three-dimensional case 
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23| . We also assume that the chosen parameters for 
the condensate are such the nonlinear energy is smaller than the zero point energy implying 
that (Ntts/az) < {Ro/azY, where Rq is the size of the cloud in the x-y plane. Throughout the 
studies we have restricted N and so that the above conditions are satisfied. Under these 
conditions we can write the condensate wave function ^ in a factorized form as following 



We take the wave function along z-axis as 

1 



(6) 



(7) 



such that 



' dz = 1. 



In order to study the time evolution of a Bose condensate confined in a Q2D trap as described 
above, we choose the general trial wavefunction in x-y coordinates as. 



(8) 



where Q;(t), and 7(t) are time dependent complex variational parameters. We have 
taken the trap to be anisotropic in the x-y plane with anisotropy parameter A given by 
ujx = \uJy = ujq\/\. In equation (jHJ, X{= x/ao) and Y{= y/ao) are dimensionless coordinate 
variables with ao = {h/mcuoY^'^. We can write the condensate parameters a(t), P{t) and 
7(t) such that 

a{t) = ai{t) + ia2{t), /3{t) = /3i{t) + i(32{t), and 7(t) = 71 (t) + i72(t). 

Here ai{t) and l3i{t) represent the inverse square of the widths along x and y-axes respec- 
tively. The imaginary parts 02 (t) and /?2(^) account for the time dependent phases. The 
parameters 71 (t) and 72 (t) representing the cross term in the coordinates account for the 



scissors modes of collective oscillation 



25l |. From the normalisation condition we can 



find that \C{t)\= NVD/n, where D = a^pi - 7^. 



By substituting equations (0)-® in equation (|T)), we can evaluate Lagrangian as 



(9) 



where 



(ai/52 + ^"2 - 27172) 



E 



an 



K 



Nhwo^ [D{ai + A) + a,{(3i + 7^) + (3,{al + 7^) - 27i72(a2 + /52)] + ^ 



E^ 



ho 



AD 



E. 



(1) 

int 



4 ^/n \a. 



+ 



= Nhw^Cial. {2na 



3/2 



^3/2^3/4_ 



To know the static properties of the condensate, we substitute the static values of different 
parameters as ai = aio, Pi = Pw, 02 = /52 = 0, 71 = 72 = in the Lagrangian in equation 
©• Then, from the Lagrange's equations of motion, we can obtain the following coupled 
equations for a 10 and /5io 

A 



22 , "-1 3/2 , <3 , ^3/4 7/4 

ao"lO - — + — Ao «10 + -k2pio "10 



'^0/^10 



^1 ,31/2 3/2 , 3 

yAo «io + 

^1 1/2^3/2 , 3 
y«io /?io + 



0, 
= 0, 



(10) 
(11) 



where 

h 

k2 



4. 



271 

128 
15V? 



an f27ra,) 



3/2 



These equations are coupled nonlinear equations. The terms with k2 as coefficient in the 
above expressions appear due to the incorporation of LHY correction in the energy e{p). 
For the known values of parameters ( N,as,az, ujq and A ) for the condensate in a Q2D trap, 
the static size-parameters aio and Piq along x and y-axes can be evaluated by solving the 
above equations. We have solved these equations numerically. The values of aio and Piq 
obtained thus have been used for the calculation of collective mode frequencies as described 
in the following section. 



III. COLLECTIVE MODES 



The low energy excitations of condensate correspond to the small oscillations in the 
cloud around the equilibrium configuration. Therefore we can expand the time dependent 
variational parameters around their static values in the following manner, 



ai = aio + 6ai{t), = + 6 (3i{t) , -f^ = Sji{t) (710 = 0), 



(12) 



After substituting the above values in the Lagrangian in equation © and using the Euler- 
lagrange equations, the time evolution of the parameters 5 Pi and ^71 can be obtained 

as 



5ai + ijJ^a^ 



- ^aioA/«io/3io - \k2aw{aiQPiof/'^ 



SPi + UqUq 



^qQq «io 
2 Pio 
4 



hiaioPioY^' + lk2iaioPiof' 



Xal 2 



6ai 
SPi = 0, 
SPi 



(13) 



^0^0 Pio 



2 aio 
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hiaiopioY''^ + -k2{aioPio] 



3/4 



571 + 



A + ^ ) + 2aQaioPio 



Sai = 0, 
571 = 0. 



(14) 
(15) 



In the above equations (fTSj) - (fT^. the first two equations (i.e. (fT^ and (IHj)) are the coupled 
equations for 6ai and 6Pi . These equations represent the monopole and quadrupole modes 
of collective oscillations. The equation (fTHjl for (57iis not coupled to the other two and 



represents the scissors mode of collective oscillation [10|. In order to calculate the the mode 
frequencies, we take the harmonic time dependence of amplitudes as 



6ai = 6ai{uj) exp{—iut), 5 Pi = SPi{uj) exp(— icut), 

571 = 6'yi{uj) exp{—iut) . 



(16) 



This results in the solution for mode frequencies as 



2 2 



ai 



+ a2 hi 



1/2 



(17) 



6 



and 



where, 



A + - ) +2a^(aio/3io), 



(18) 



bi 



2 /3: 



10 



fci(aio/3io)'/' + ^A;2(aio/3io)'/" 



02 



2 aio 



9 

/ci(aioAo)^'^^ + |/j2(aio/3io! 



3/4 



62 = cooal 



ki 



- yAoVW^- g/i;2/3io(Qiio/3io)^^^ 
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FIG. 1: Variation of monopole mode frequency with scattering length. The other parameters are 
N = 10^, a;^ = 27r X lOOOHz, loq = 27: x 2QHz, a„ = 27.5^° and A = 0.8. The dotted curve shows 
frequency values calculated within the mean-field theory (i.e. k2 = 0) whereas the continuous curve 
shows that after incorporating LHY correction term (i.e. ^2 = finite). 

Here a;+ is monopole mode frequency, a;_ is quadrupole mode frequency and cUs is scissors 
mode frequency. The terms with k2 as coefficient in the above expressions of frequencies 
represent the corrections due to LHY term in the interaction energy £{p). It can be noted 
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that monopole and quadrupole mode frequencies depend explicitly upon k2 whereas scissors 
mode frequency does not depend upon it. Variation of monopole mode frequency uj+ with 
scattering length (i.e. ag/cLn, a„ is scattering length in absence of magnetic field) is shown 
in FIG. n It is evident that with finite value of ^2, monopole mode frequency changes with 
as/ an in a differnt manner than that with k2 = 0. The frequency values from the two curves 
differ in quantitative as well as in qualitative manner as shown in the figure. After the initial 
fall in the frequency uj^ with increasing ag/an, it becomes nearly constant if ^2 = 0. However 
for a finite value of k2, after the initial fall, there is monotonous increase in values with 
increasing ag/an- Although the modifications in the values of 0;+ with finite ^2 are ius t b\ 
a few percent, nonetheless, they can be easily detected as has been pointed out earlier 
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FIG. 2: Variation of mode splitting ( A = ujs — 0J~ ) with as/a-n- The other parameters are 
= 27r X lOOOiJz, ojq = 2tt X 20Hz, V = 10^, a„ = 27.5A° and A = 0.8. The dotted curve shows 
A values calculated within the mean-field theory (i.e. /c2 = 0) whereas the continuous curve shows 
that with incorporating LHY term in the inteaction energy (i.e. k2 = finite). 

As has been discussed earlier [l^ the splitting between quadrupole and scissors modes 
occurs in a 2D anisotropic trapped Bose condensate. The value of mode splitting (A = 
— uj^) depends upon anisotropy and interaction strength. To assess the effect of LHY 
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term on the value of A, we have calculated it as function of as /an for finite ^2 as well as for 
k2 = 0. The results are as shown in the FIG. EJ The values shown here are for anisotropy 
parameter A = 0.8. As is evident from the figure, the splitting A calculated with finite k2 
keeps on changing with as/ an whereas it becomes nearly constant with as/ an for k2 = 0. It 
is worth mentioning that the scissors mode frequency Us does not explicitly depend upon 
the LHY term ( i.e. ^2)- However, its dependence on this term occurs via parameters aio 
and jSio, which is quite small. Therefore the difference in the A values in two curves in FIG|21 
arises mainly due to dependence of quadrupole frequency uj_ on the terms containing /c2- 

IV. CONCLUSION 

We have theoretically studied the collective mode frequencies of a quasi two-dimensional 
Bose condensate in the large gas parameter regime. We find that in this regime, the in- 
clusion of LHY correction term beyond the mean-field term in the interaction energy leads 
to quantitative as well as qualitative changes in the dependence of collective mode frequen- 
cies on interaction strength. Such regime of parameters is experimentally accessible with 
the present advancements in this field and modifications in the mode frequencies are easily 
detectable. 
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